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Abstract 

A shape phase transition is demonstrated to occur in 190 W by applying the Pro- 
jected Shell Model, which goes beyond the usual mean-field approximation. Rotation 
alignment of neutrons in the high-j, ? 13 / 2 orbital drives the yrast sequence of the 
system, changing suddenly from prolate to oblate shape at angular momentum IQh. 
We propose observables to test the picture. 
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Understanding the fundamental excitations of many-fermion systems and tran- 
sitions between different excitation modes is an important issue in mesoscale 
physics. In nuclei where the ground state is typically a superconducting state 
with paired nucleons moving in orbits, the low-lying excitation spectrum is 
generally formed by collective motion (for example, rotation and vibration 
with different nuclear shapes, or deformation) and nucleon pair breaking. Be- 
ing a finite-size quantum system, for certain numbers of protons and neutrons, 
a subtle rearrangement of only a few nucleons among the orbitals near the 
Fermi surface can result in completely different collective modes, constituting 
a suitable situation for studying the shape phase transition [Tf2"f5] . 

There is evidence for several different types of shape phase transition in nuclei. 
Casten et al. [I] systematically studied spherical-to-deformed ground-state 
phase transitions as a function of neutron number and proton number; and 
Regan et al. [5] observed the evolution from collective vibrational to rotational 
structure as a function of angular momentum (see also theoretical discussions 
[6117] ). What we discuss in this Letter is a robust example of a shape phase 
transition that occurs along an yrast cascade (the locus of lowest energy states 
of each spin) between states of prolate and oblate shape in an isolated nucleus. 

The neutron- rich Hf and W nuclei with neutron number iV »s 116 are expected 
to exhibit interesting critical-point phenomena due to competing collectivity 
of prolate and oblate shapes. Since experimental evidence is sparse, we refer 
mainly to theoretical predictions (e.g. Refs. [SIS])- In these nuclei, the ground 
state and low-spin states are of a prolate shape. However, the occupation of 
low-f^, i\3/2 neutron intruder orbits when the shape is oblate favors rotation 
alignment at moderate spin values. (K is the symmetry-axis projection of the 
angular- momentum vector.) This was suggested by energy-surface calculations 
for 190 W based on a mean- field method [TU]. Ref. fTDj also suggested that the 
190 W structure is a candidate for a classic case of "giant backbending" , of the 
type originally described by Hilton and Mang [TT] in the 1970's. Confirma- 
tion of these suggestions requires additional theoretical investigations as well 
as experimental data. Here we show, with aid of detailed shell-model-type 
calculations, that 190 W can exhibit an excellent example of a shape phase 
transition that differs from the previously known types of phase transition 
reported in the literature [Ifo] . The combination of electrical quadrupole mo- 
ments or B(E2) values, and gyromagnetic ratios (g factors), should be able to 
give firm confirmation. The isomerism [T21T3] is a key feature in giving access 
to these observables. 

Phase transitions in nuclei have been theoretically described mainly by means 
of algebraic models [6fT][H] . in which a transition-driving control parameter 
appears explicitly in the Hamiltonian. There have been microscopic studies 
based on mean-field theories (for example, Ref. [15]). Using the Landau the- 
ory, Alhassid et al. [IS] found a universal behavior of rapid changes of the 
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Fig. 1. Energy surface calculation for 190 W. (a) Energies obtained with projection 
on qp vacuum states only, (b) Energies from calculation with basis that includes 
2-qp states. 

equilibrium shape from the prolate ground state to oblate excited states in 
hot rotating nuclei. However, the nuclear shell model can provide a more fun- 
damental basis to study phase transitions. The advantage of a shell-model 
study is that one may see the microscopic origin of a phase transition by ana- 
lyzing wave functions. The early work of Federman et al. [17 \ showed perhaps 
the first example of such a kind. The more recent phase transition study with 
the Monte Carlo Shell Model [18] and large-scale shell-model calculations [19] 
are other examples. However, these shell models are applicable only to small 
nuclear systems, and 190 W is far too large for meaningful calculations of this 
type. 

The Projected Shell Model (PSM) [20] is a shell model that uses deformed 
bases and the projection technique. It is applicable to any large size of de- 
formed systems (with extreme examples from superdeformed [21] to super- 
heavy nuclei [22J). The PSM's two-body residual interactions are of the quadrupole 
plus pairing type, with the quadrupole-pairing term included |20j. Wave func- 
tions of the PSM are written in terms of angular-momentum-projected multi- 
quasiparticle (qp) states 

iM-T.fAM, (i) 

K 



where the index k labels basis states. P^k * s ^ ne angular-momentum-projection 
operator [20] and the coefficients are weights of the basis states. For even- 
even nuclei, \4> K ) in Eq. ([I]) is 

{|0),atat|0),44|0),---} (2) 
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where aj, and a\ are the creation operator for neutrons and protons, respec- 
tively. The qp states are obtained from a deformed Nilsson calculation (with 
the Nilsson parameters taken from Ref. [23]) followed by a BCS calculation, 
in a model space with three major shells for each kind of nucleon (N = 4, 5, 
6 for neutrons and N = 3, 4, 5 for protons). The corresponding qp vacuum is 
|0) = \e) at deformation e. We assume axial symmetry in the basis (consistent 
with Ref. [10]) and, therefore, each basis state in ([2]) can be labeled by the K 
quantum number. While the basis states have good K, the states of Eq. flTJ) 
generally do not, as they are linear combinations of various if-states. 

Fig. 1 shows the projected energy surface calculation for 190 W. The quantities 
plotted in Fig. 1 are 



These are energies with different angular momenta I — 0, 2, • • ■ calculated as a 
function of basis deformation e, varying from negative values (corresponding to 
oblate shapes) to positive values (corresponding to prolate shapes). The energy 
curves show that there are two pronounced minima, sitting respectively at the 
prolate and oblate side, with |e| m 0.15. The minimum at e ~ +0.15 is lower, 
suggesting that the ground state of 190 W is of prolate shape, which is consistent 
with the ground state deformation obtained in Refs. PfTO] . The excited 1 = 
state, i.e. the minimum at the oblate side with e ~ —0.15, is about 0.6 MeV 
higher than the ground state. Thus the result predicts a competing picture 
of prolate and oblate shapes with two sets of corresponding collective states. 
Dynamic perturbations may change the balance in the competition as the 
order parameter varies, causing a prolate-to-oblate phase transition. 

Perturbations in quantum systems can arise from different sources. The dy- 
namic process responsible for the phase transition in this case is the nuclear 
rotation and the order parameter is the total spin /. As the nucleus rotates, 
individual nucleon pairs tend to align their spins with the collective rotation 
axis through the Coriolis force. The alignment lowers the system's energy; the 
amount of energy gain depends primarily on the microscopic structure of the 
aligned particles. To see this clearly, we compare two calculations in Fig. 1. 
The left figure shows the calculation in which contains the projected qp 
vacuum state only and all 2-qp states in basis ([2]) are switched off. The under- 
lying physics is that there is no rotation alignment allowed in the calculation. 
In such a case, minima on the prolate side remain lower for every state. This 
implies that the prolate shape always wins in the competition with the oblate 
shape if no rotation alignment can happen. 

Calculations that include 2-qp states in the model space are shown in the 
right figure in Fig. 1. With inclusion of 2-qp states, energy levels of higher 
spins are compressed considerably. As the arrows indicate, along the yrast 
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Fig. 2. Selected bands in the 190 W calculation, (a) Bands with oblate shape (calcu- 
lated at e = —0.15). (b) Bands with prolate shape (calculated at e = +0.15). 

cascade, the path for the lowest state at each / starts from the prolate ground 
state 1 = 0, goes up till / = 10, and then jumps to the oblate side. The wave 
functions with prolate and oblate shapes are expected to be very different, with 
a very small overlap between them. Therefore, the above result implies that 
the yrast cascade of 7-rays is broken and the 7 transition between the oblate 
and prolate states may be strongly hindered. This is a favorable situation for 
the oblate state at the phase transition to be isomeric. We emphasize that 
this type of calculation contains beyond-mean-field correlations in the sense 
that each point in the curves in Fig. 1 is obtained by projecting the intrinsic 
configurations (deformed mean-field results) onto states with good angular 
momentum and mixing the projected states through residual interactions. 

A phase transition in a quantum many-body system generally refers to an 
abrupt, qualitative change in the wave function. This subject is of significant 
interest for many subfields, and is an attractive current topic in nuclear physics 
[T1I21I3] . Our results in Fig. 1 show a clear example of a shape phase transition 
occurring in an isolated nucleus between two deformed (prolate and oblate) 
shapes driven by the 113/2 rotation alignment. Coexistence of near-spherical 
and deformed shapes is a known effect for nuclei near the proton shell closure 
(see for example 186 Hg [24)). However, the present case with N = 116, which 
goes further to the neutron-rich region and has additional valence protons, 
is more robust. In 186 Hg the ground band is close to spherical and not well 
developed, but there is also a spin- 16 band crossing. This latter feature has 
parallels with the behavior of 190 W, where the prolate (oblate) wave function 
contains the highest (lowest) K-components of the neutron i 13 / 2 orbit (see 
later discussions), and the wave functions before and after the transition are 
therefore sharply different. 
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As long as a nucleus exhibits well-defined deformation minima, such as those 
shown in Fig. 1, one can apply the PSM and perform detailed shell-model 
analysis at the minima. To further understand the results in Fig. 1, we show 
several representative bands in Fig. 2. Angular-momentum-projection on a 
multi-qp state \4> K ) with a sequence of / generates a band. One may define the 
rotational energy of a band (band energy) using the expectation values of the 
Hamiltonian with respect to the projected \<f) K ) [20] 



These are bands before configuration mixing (diagonal elements in the Hamil- 
tonian matrix) and are denoted as "bef diag" in Fig. 2. Those after configu- 
ration mixing, i.e. the superposition of all basis states k in (JTJ), are denoted 
as "aft diag" . Bands after mixing are solutions of the eigenvalue problem and 
dynamics is thus taken into account. 

In the calculation for Fig. 2, shell model diagonalization is carried out at two 
deformations, e = ±0.15 (see Fig. 1). Note that in the calculation, about 100 
2-qp states are included in each diagonalization. We present only a few of them 
to illustrate the physics. We recall that in 190 W, the neutron occupation is very 
different at the two energy minima. On the prolate side, one can find high-.fr 
orbits such as §[624] and y[615] close to the neutron Fermi surface, while 
on the oblate side, only the lowest K-states of the zi 3 / 2 intruder shell appear 
around the Fermi surface. In the right part of Fig. 2, the ground band (dotted 
curve) starts from the origin and goes up monotonically. At I = 6, the K = 1, 
2-qp band (dotted-dashed) built by | [624] © y[615] quasi- neutrons crosses the 
ground band. However, interactions at the crossing region are so strong that 
they repel the crossing bands from each other, forming two smooth sequences. 
In fact, the resulting bands (solid and dashed curves) are smooth bands as if 
no band crossing had occurred. The lowest states (solid curve) obtained after 
configuration mixing reproduce very well the known data (filled diamonds) 
jl3j . On the other hand, in the left part of Fig. 2, the behavior of the low-K 
2-qp states on the oblate side, §[660] ©§[651] and §[660] ©§[642] (both dotted- 
dashed), is quite different from the 2-qp state on the prolate side. After 1 = 6, 
the combined behavior of these two bands gives rise to nearly degenerate states 
at / = 8, 10, and 12 lying lowest in energy. These configurations are predicted 
to be dominant in the lowest oblate states (solid curve), which compares well 
with the tentative data point from Ref. [13]. 

One can thus understand that the compression of the states with 1 = 8, 10, and 
12 of the oblate shape is a consequence of the particular rotation-alignment 
effect of the \ow-K, i 13 / 2 neutrons. The effect considerably lowers the total 
energy, which helps the oblate states to eventually win the prolate-oblate 
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Fig. 3. Predicted electromagnetic properties for 190 W. (a) B(E2, I — > I — 2) values, 
(b) g factors. 

competition. As a result, a prolate-to-oblate shape phase transition occurs 
along the yrast cascade, and the oblate state with I = 10 becomes isomeric. 

Though not shown in Fig. 2, we have found several 2-qp high-iT bands with 
prolate shape in the calculation. The bandhead energies are ranging down 
to 1.2 MeV of excitation. As some of them lie low in energy and have K- 
quantum numbers much differing from the K = ground band, these high-.fr 
bands are also predicted to be X-isomers. The dominant structures of the 
lowest such bands are: K w = 6+ §[512] © §[505]; K n = 10+, §[624] © y[615]; 
K n = 7", §[512] © y [615]; and K* = 10", §[505] © y[615]. 

The energy minimum with oblate shape (see Fig. 1), which lies about 0.6 MeV 
above the ground state, is a + shape isomer predicted in our calculation. A 
+ shape isomer has been observed in 72 Kr [23]. A recent experiment at GSI 
may have found evidence [26] for /3-decay to this shape isomer. 

A crucial test for the above predictions is to measure additional observables 
for the states before and after the phase transition. As the states belong to 
very different structures, one may find clear indications when studying their 
electromagnetic properties, e.g. B(E2) values and g-factors. B(E2) values that 
are related to the electric quadrupole transition probability from an initial 
state 7 to a final state 1 — 2 are given by 

B(E2, 1 - / - 2) = ^ | ( 11411^ ) I 2 • (5) 



7 



The effective charges used in the calculation are the standard ones: e n = 1.5e 
and e v = 0.5e. Magnetic properties are described through g factors, which are 
defined as 

g{1) = M (6) 

fi N Jl(I + l) 



In the (/-factor calculation, we use the standard values for gi and g s , namely, 
the free values for gi and the free values damped by a 0.75 factor for g s . In 
Eqs. (IS]) and (jH]), wave functions ip 1 ) are those of Eq. ([T]). 



Presented in the upper part of Fig. 3 are two sets of B(E2) values, calculated 
at the prolate and oblate minima. It can be seen that the values at the oblate 
minimum are clearly predicted to be smaller than those at the prolate min- 
imum for the spin states / < 6. A drastic drop of the oblate B(E2) occurs 
at / = 8, due to the sharp band crossing of the low-K z 13 / 2 neutrons. Thus 
we predict very different B(E2) values (and hence also quadrupole moments) 
for the prolate and oblate structures. Similar conclusions can be drawn for 
g factors. As the lower part of Fig. 3 shows, there are significant differences 
between the prolate and oblate g factors. They both show a drop at / = 6 
and 8, due to the rotation alignment of neutrons; however a much greater 
drop occurs for the oblate g factors, which even causes the sign to change. 
Experimental observation of the differences in electromagnetic properties will 
definitely determine the structure around the crossing point, and thus provide 
a crucial test for our prediction. We note that new experimental information 
may soon be available 



In conclusion, by using the Projected Shell Model, which conserves total angu- 
lar momentum and includes configuration mixing beyond the usual mean-field, 
we have demonstrated a new type of shape phase transition that occurs along 
the yrast cascade in 190 W. The transition is driven by the dynamic process 
of rotation alignment of the high-j neutrons, which leads to the preferred 
state suddenly changing shape from prolate to oblate at / = 10. The tran- 
sition results in oblate shape isomers, extending the predictions of Ref. [TO] , 
We have proposed testable quantities for future experiments to measure. Fi- 
nally, we note that our searching for energy minima has been restricted to 
axial symmetry. A more general case will be minimization in the e-7 plane 
with three-dimensional angular momentum projection. This is the subject of 
future work. 
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